Reaction rates for the s-process neutron source 22 Ne+a 
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The 22 Ne(a,n) 25 Mg reaction is an important source of neutrons for the s-process. In massive stars responsi- 
ble for the weak component of the s-process, 22 Ne(a,n) 25 Mg is the dominant source of neutrons, both during 
core helium burning and in shell carbon burning. For the main s-process component produced in Asymptotic 
Giant Branch (AGB) stars, the 13 C(a,n) 16 reaction is the dominant source of neutrons operating during the 
interpulse period, with the 22 Ne+a source affecting mainly the s-process branchings during a thermal pulse. 
Rate uncertainties in the competing 22 Ne(a,n) 25 Mg and 22 Ne(a,7) 26 Mg reactions result in large variations of 
s-process nucleosynthesis. Here, we present up-to-date and statistically rigorous 22 Ne+a reaction rates using 
recent experimental results and Monte Carlo sampling. Our new rates are used in post-processing nucleosyn- 
thesis calculations both for massive stars and AGB stars. We demonstrate that the nucleosynthesis uncertainties 
arising from the new rates are dramatically reduced in comparison to previously published results, but several 
ambiguities in the present data must still be addressed. Recommendations for further study to resolve these 
issues are provided. 

PACS numbers: 26.20.Kn - 26.20.Fj - 25.55.-e - 24.30.-v 



I. INTRODUCTION 

The s-process is responsible for creating about half of the 
elements heavier than iron that are observed in the solar sys- 
tem OJJ] . This process involves the slow capture of neutrons 
(slower than the average /3-decay rate of unstable nuclei) onto 
seed material, hence nucleosynthesis follows the nuclear val- 
ley of stability. By considering the solar system abundances 
of s-only nuclei (that is, nuclei that can only be produced in 
the s-process) it can be shown that there are two key compo- 
nents of the s-process: the "main" component and the "weak" 
component J2|]. The main component produces s-nuclei with 
masses of A > 90, while the weak component enriches the 
s-nuclei abundances at A < 90. 

The main component of the s-process arises from neutron 
captures during He-buming in M < AMq Asymptotic Giant 
Branch (AGB) stars (a detailed discussion of nuclear burn- 
ing in AGB stars can be found in Refs. d and 01). In 
low mass (0.8 to 4 M Q ) AGB stars of solar metalicity, most 
neutrons are released through the 13 C(a,n) 16 reaction dur- 
ing the inter-pulse period, while the 22 Ne(a,n) 25 Mg reac- 
tion produces an additional burst of neutrons during thermal 
pulses. This burst of neutrons affects mainly the branchings 
in the s-process path. In intermediate-mass AGB stars (M > 
4Mq), where the temperatures are expected to be higher, the 
22 Ne(a,n) 25 Mg reaction is thought to be the main source of 
neutrons and could explain the enhancement of rubidium seen 
in some metal poor AGB stars iH-IH]. In addition to s-process 
elements, the 22 Ne(a,n) 25 Mg and 22 Ne(a,7) 26 Mg rates influ- 
ence the relative production of 25 Mg and 26 Mg, whose abun- 
dance ratios can be measured to high precision in circumstel- 
lar ("presolar") dust grains. Magnesium is also one of the 
few elements for which the isotopic ratios ( 25 Mg/ 24 Mg and 
26 Mg/ 24 Mg can be derived from stellar spectra (for example, 
Refs. ITolflill ). However, Karakas et al. Il2tl showed that with 



their estimated 22 Ne(a,n) 25 Mg and 22 Ne(a,7) 26 Mg reaction 
rate uncertainties, the relative abundances of 25 Mg and 26 Mg 
predicted by their stellar models can vary by up to 60%. 

The weak component of the s-process arises from nuclear 
burning in massive stars. The core temperature in these stars 
(typically with M > 11M Q ) becomes high enough during 
He-burning for the 22 Ne(a,n) 25 Mg reaction to produce a high 
flux of neutrons shortly before the helium fuel is exhausted. 
Any remaining 22 Ne releases a second flux of neutrons dur- 
ing convective carbon shell burning. The s-process yield in 
these stars is therefore sensitive to the temperature at which 
the 22 Ne(a,n) 25 Mg reaction starts to produce an appreciable 
flux of neutrons. The et al. IU3I1 showed that the s-process 
during the core He-buming stage in massive stars depends 
strongly on both the 22 Ne+a and the 16 0(n,7) 17 reaction 
rates. They also found that not only are the overall uncer- 
tainties in the rates important, but also the temperature depen- 
dence of the rates. 

The 22 Ne+a reactions also affect nucleosynthesis in other 
astrophysical environments. During type II supernova explo- 
sions, two 7-ray emitting radionuclides, 26 A1 and 60 Fe are 
ejected, and their abundance ratio provides a sensitive con- 
straint on stellar models [14, and references therein]. The 
species 60 Fe is mainly produced in massive stars by neutron 
captures during convective shell carbon burning [e.g.,[l5l. Its 
abundance, therefore, depends strongly on the 22 Ne+a rates. 
The 22 Ne+a rates also play a role in type la supernovae. 
Throughout the "simmering" stage, roughly 1000 years prior 
to the explosion, Piro and Bildsten [16] suggested that neu- 
trons released by the 22 Ne(a,n) 25 Mg reaction affect the car- 
bon abundance, thus altering the amount of 56 Ni produced 
(i.e. , the peak luminosity) in the explosion. Timmes et al. 
IU7I1 also found that during the explosion, neutronisation by 
the 22 Ne(a,n) 25 Mg reaction affects the electron mole fraction, 
Y e , thus influencing the nature of the explosion. 
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In this work we will evaluate new reaction rates for 
22 Ne+a . Compared to previous results 1IT2I [l8l [l9ll our new 
rates are significantly improved because (i) we incorporate all 
the recently obtained data on resonance fluorescence absorp- 
tion, a-particle transfer etc., and (ii) we employ a sophisti- 
cated (Monte Carlo) method to estimate the rates and associ- 
ated uncertainties. We have recently presented new 22 Ne+a 
rates in Ref. [20], but did not give a detailed account of their 
calculation. Since the latter results were published, we found, 
and could account for, a number of inconsistencies in data pre- 
viously reported in the literature. In addition, new data from 
Ref. Il2lll became available, which have been included in the 
present work. Thus the rates presented here supersede our 
earlier results 1I20I1 . 

The paper will be organised as follows: in Sec.[n]a detailed 
discussion of the Monte Carlo method used to calculate re- 
action rates is discussed. This method is described in detail 
elsewhere I22I1 but will be summarised to show its applicabil- 
ity to the specific cases of the 22 Ne+a reactions. The 22 Ne+a 
rate calculations and comparisons with the literature will be 
presented in Sec. [HI] The reaction rates will then be used to 
present new nucleosynthesis yields along with their uncertain- 
ties in Sec.[VT] Conclusions will be presented in Sec. I VII I 



II. REACTION RATE FORMALISM 
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where the resonance strength, uij is defined by 

w) = LoT a Tb/T, 

E r is the resonance energy; T a (E), Tb(E), and T(E) are the 
energy-dependent entrance channel (particle) partial width, 
exit channel partial width, and total width, respectively; and 
uj, the statistical spin factor, is defined by to = (2J + 
l)/(2Jo + l)(2Ji + 1), where J and J,; are the resonance 
and particle spins, respectively. The particle partial width, T c , 
can be written as the product of an energy-independent re- 
duced width, 7 2 , and an energy-dependent penetration factor, 
P c (E),as 



r c = 2P C (£) 7 2 . 



(5) 



For the present case of 22 Ne+a, the entrance channel (a- 
particle) reduced width, 7 2 , is related to the a-particle spec- 
troscopic factor, S a , by 



A. Thermonuclear Reaction Rates 

The reaction rate per particle pair in a plasma of tempera- 
ture, T, is given by 
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where /1 is the reduced mass of the reacting particles, /i = 
MqM\/{Mq + Mi)', Mi denotes the masses of the particles; 
k is the Boltzmann constant; E is the centre-of-mass energy 
of the reacting particles; and <r(E) is the reaction cross section 
at energy, E. 

The strategy for determining reaction rates from Eq. ((T) de- 
pends on the nature of the cross section. In many cases the 
cross section can be separated into non-resonant and reso- 
nant parts. Reactions such as 22 Ne+a proceed through the 
compound nucleus 26 Mg at relatively high excitation energy 
(Q Q7 = 10614.787(33) keV JH) and are frequently domi- 
nated by resonant capture. The non-resonant part of the cross 
section will, therefore, be neglected in the following discus- 
sion. The reader is referred to Refs. [l24ll and ll22l for more 
details. 

The resonant part of the cross-section can be represented 
in one of two ways: (i) by narrow resonances, whose par- 
tial widths can be assumed to be approximately constant over 
the resonance width ("narrow resonances"), and (ii) by wide 
resonances, for which the resonant cross section must be in- 
tegrated numerically to account for the energy dependence of 
the partial widths involved. The reaction rates per particle pair 
for single, isolated narrow and wide resonances, respectively, 
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where 0(a) is the single-particle radial wave function at the 
channel radius, a [see, for example, The constant 

H 2 /(/j,a 2 ) is the Wigner Limit (in the notation of Lane and 
Thomas 112711 ). It can be regarded as an upper limit, according 
to the sum rules in the dispersion theory of nuclear reactions, 
i-e., # 2 < 1. Note that it is frequently assumed that 2 = S, 
which must be regarded as a crude approximation only. For 
example, in the case of 17 levels it was shown in Ref. l28tl 
that S a exceeds B 2 a by at most a factor of 2. We will return to 
these issues in Sec. IIIICI 

The above relationships are useful since they allow for 
an estimation of the important a-particle partial widths from 
spectroscopic factors obtained in a-particle transfer reactions, 
as will be discussed later. It is important to note that the value 
of S a depends on the parameters of the nuclear potentials as- 
sumed in the transfer data analysis. Similarly, the value of 7^ 
depends on the channel radius. However, if, throughout the 
analysis, consistent values of these parameters (such as a) are 
used, their impact on the value of r Q will be strongly reduced. 



B. Monte Carlo Reaction Rates 



The equations outlined in Sec. Ill Al provide the tools for 
calculating thermonuclear reaction rates given available esti- 
mates for the cross section parameters (E r , ujj, etc.). A prob- 
lem arises, however, when statistically rigorous uncertainties 
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of the reaction rates are desired. What is usually presented in 
the literature are recommended rates, together with upper and 
lower "limits", but the reported values are not derived from a 
suitable probability density function. Therefore, the reported 
values have no rigorous statistical meaning. An attempt to 
construct a method for analytical uncertainty propagation of 
reaction rates was made by Thompson and Iliadis [29|. How- 
ever, their method is applicable only in special cases, when 
the uncertainties in resonance parameters are relatively small. 
Thompson and Iliadis 12911 were also not able to treat the un- 
certainty propagation for reaction rates that need to be inte- 
grated or for rates that include upper limits on some param- 
eters. For these reasons, a Monte Carlo method is used in 
the present study to calculate statistically meaningful reaction 
rates. 

The general strategy of Monte Carlo uncertainty^ propa- 
gation is as follows: (i) randomly sample from the probabil- 
ity density distribution of each input parameter; (ii) calculate 
the reaction rates for each randomly sampled parameter set 
on a grid of temperatures (using the same set at each temper- 
ature); (iii) repeat steps (i)-(ii) many times (on the order of 
5000). Steps (i)-(iii) will result in a distribution at each tem- 
perature grid point that can be interpreted as the probability 
density function of the reaction rate. Extraction of uncertain- 
ties from this distribution will be discussed later. While input 
parameter sampling is being performed, care must be taken 
to consider correlations in parameters. For example, parti- 
cle partial widths depend on the penetration factor, which is 
an energy dependent quantity. The individual energy samples 
must, therefore, be propagated consistently through resonance 
energy and partial width estimation in order to fully account 
for the correlation of these quantities. The code RatesMC 
l22tl was used to perform the Monte Carlo sampling and to 
analyse the probability densities of the total reaction rates. 

In order to apply Monte Carlo sampling to calculate reac- 
tion rate uncertainties, sampling distributions must be chosen 
for each input parameter. Once a reaction rate (output) dis- 
tribution has been computed, an appropriate mathematical de- 
scription must be found to present the result in a convenient 
manner. Statistical distributions important for reaction rate 
calculations are described in detail in Refs. Il22ll and 130J, and 
are summarised briefly below. 

Uncertainties of resonance energies are determined by the 
sum of different contributions. In this case, the central 
limit theorem of statistics predicts that resonance energies are 
Gaussian distributed. Note that there is a finite probability of 
calculating a negative resonance energy and that this choice 
of probability density naturally accounts for the inclusion of 
sub-threshold resonances in the above formalism. A reso- 
nance strength or a partial width, on the other hand, is exper- 
imentally derived from the product of measured input quanti- 
ties (e.g., count rates, stopping powers, detection efficiencies, 



1 Throughout this work, care is taken to refer to the terms uncertainty and 
error correctly. The term error refers to a quantity that is believed to be 
incorrect, whereas uncertainty refers to the spread estimate of a parameter. 



etc.). In such a case the central limit theorem predicts that res- 
onance strengths or partial widths are lognormally distributed. 

The lognormal probability density for a resonance strength 
or a partial width is given by 
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with the lognormal parameters /i and a representing the mean 
and standard deviation of In x. These quantities are related to 
the expectation value, E[x], and variance, V[x], by 
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The quantities, ln(2?[x]) and can be associated with 

the central value and uncertainty, respectively, that are com- 
monly reported. Note that a lognormal distribution is only 
defined for positive values of x. This feature is crucial be- 
cause it removes the finite probability of sampling unphysical, 
negative values when Gaussian uncertainties are used. This 
is especially true for partial width measurements, which fre- 
quently have uncertainties in the 20-50% range. Note, also, 
that a 50% Gaussian uncertainty results in a 3% probability of 
the partial width having a value below zero. 

The important problem of estimating reaction rates when 
only upper limits of resonance strengths or partial widths are 
available will now be discussed. The standard practise in nu- 
clear astrophysics [see, for example. Il8[ l3lTl is to adopt 10% 
resonance strength upper limit values for the calculation of the 
recommended total rates. "Lower limits" or "upper limits" of 
rates are then derived by completely excluding or by adopting 
the full upper limit, respectively, for all resonance strengths. 
This procedure is questionable for two reasons. First, without 
further knowledge, it is implicitly assumed that the probabil- 
ity density for the resonance strength is a uniform distribution 
extending from zero to the upper limit value. The implication 
is that the mean value of the resonance strength amounts to 
half of its upper limit value. This conclusion contradicts fun- 
damental nuclear physics, as will be explained below. Sec- 
ond, the derived "upper limit" and "lower limit" on the total 
reaction rate are usually interpreted as sharp boundaries. This 
conclusion is also unphysical, as will be explained below. 

The strength of a resonance depends on particle par- 
tial widths, which can be expressed in terms of reduced 
widths, 7 2 , or, alternatively, spectroscopic factors, S (see sec- 
tion III At . These quantities depend on the overlap between 
the incoming channel (a + A) and the compound nucleus 
final state, which in turn depends on a nuclear matrix ele- 
ment. If the nuclear matrix element has contributions from 
many different parts of configuration space, and if the signs of 
these contributions are random, then the central limit theorem 
predicts that the probability density of the transition ampli- 
tude will tend toward a Gaussian distribution centred at zero. 
The probability density of the reduced width, representing the 
square of the amplitude, is then given by a chi-squared dis- 
tribution with one degree of freedom. These arguments were 
first presented by Porter and Thomas 113 211 and this probability 
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density is also known as the Porter-Thomas distribution. For 
a particle channel it can be written as 



3 -e 2 /(20 ) 



(10) 



where c is a normalisation constant, 8 2 is the dimensionless re- 
duced width, and 8 2 is the local mean value of the dimension- 
less reduced width. The distribution implies that the reduced 
width for a given nucleus and set of quantum numbers varies 
by several orders of magnitude, with a higher probability the 
smaller the value of the reduced width. The Porter- Thomas 
distribution emerges naturally from the Gaussian orthogonal 
ensemble of random matrix theory and is well established ex- 
perimentally (see Ref. ll33ll for a recent reviewfl 

The above discussion provides a physically sound method 
for randomly sampling reduced widths (or spectroscopic fac- 
tors) if only an upper limit value is available. Furthermore, in 
the present work we assume a sharp truncation of the Porter- 
Thomas distribution at the upper limit value for the dimen- 
sionless reduced width, 9 2 d , that is, we randomly sample over 
the probability density 
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Once dimensionless reduced widths are obtained from sam- 
pling according to equation (fTTT i. samples of particle partial 
widths can be found from equation (O. Subsequently, samples 
of resonance strengths can be determined from equation (|4}. 

In order to utilise equation (fTTT i for Monte Carlo sampling 
of a-particle partial widths, the mean value of the dimension- 
less reduced width, 8 2 a , must be known. To this end we con- 
sidered 360 a-particle reduced widths in the A=20-40 mass 
region [see 36, and references therein]. The distribution is 
shown in figure Q] as a black histogram. Binning and fit- 
ting the data to equation ( TTOb (solid line) results in a best-fit 
value of 8 2 a = 0.010, which we adopt in the present work. 
It is important to recall the above arguments: the distribution 
of reduced widths for a given nucleus, given orbital angular 
momentum, given channel spin, etc., is expected to follow 
a Porter- Thomas distribution. However, because of the rel- 
atively small sample size of 360 values, we were compelled 
to fit the entire set by disregarding differences in nuclear mass 
number and orbital angular momentum. For this reason, our 
derived mean value of 0.010 must be regarded as preliminary. 
More reliable estimates of 8 2 a have to await the analysis of a 



2 Recently, a high precision study of neutron partial widths in plutonium 
(A = 192, 194) by Koehler et al. and a re-analysis of the Nuclear 
Data Ensemble (A = 64 — 238) in Ref. f35ll have claimed that the data are 
not well described by a \ 2 distribution with one degree of freedom (u = 1, 
i.e., a Porter-Thomas distribution). They find, depending on the data set 
under consideration, values between v = 0.5 and v = 1.2. These new 
results are controversial and more studies are needed before the issue can 
be settled. It is not clear at present if this controversy has any implications 
for the compound nucleus 26 Mg. 



significantly larger data set of a-particle reduced widths when 
it becomes available in the future. 

From the arguments presented above it should also be clear 
that the Porter-Thomas distribution is not expected to repre- 
sent the reduced width of all nuclear levels, particularly if the 
amplitude is dominated by a few large contributions of config- 
uration space. The most important example for the latter situa- 
tion are a-cluster states, which are expected to have relatively 
large reduced widths. Indeed, the large reduced width values 
in figure[T]that are not described by the Porter-Thomas distri- 
bution (solid line) originate most likely from a-cluster states. 
Clearly, the nuclear structure of a level in question must be 
considered carefully. For this reason, results from a-particle 
transfer studies are very important. It can be argued that these 
measurements populate preferentially a-cluster states, with 
large reduced widths (or spectroscopic factors), while levels 
not populated in a-transfer have small reduced widths and, 
therefore, are more likely statistical in nature (i.e., described 
by a Porter-Thomas distribution). This issue will become im- 
portant in later sections. 

Once a random sampling of all input parameters has been 
performed, an ensemble of reaction rates is obtained. From 
its probability density one can extract descriptive statistics 
(mean, median, variance etc.). For the recommended reaction 
rate, we adopt the median value. The median is a useful statis- 
tic because exactly half of the calculated rates lie above this 
value and half below. Note that we do not use the mean value 
because it is strongly affected by outliers in the reaction rate 
distribution. The low and high reaction rates are obtained by 
assuming a 68% coverage probability. There are several meth- 
ods for obtaining these coverage probabilities, such as finding 
the coverage that minimises the range of the uncertainties, or 
one that is centred on the median. In the present work, the 16 th 
to 84 th percentiles of the cumulative reaction rate distribution 
are used. We emphasise an important point regarding reaction 
rate uncertainties: contrary to previous work, our "low" and 
"high" rates do not represent sharp boundaries (i.e., a prob- 
ability density of zero outside the boundaries). As with any 
other continuous probability density function, these values de- 
pend on the assumed coverage probability, i.e., assuming a 
larger coverage will result in a larger uncertainty of the total 
reaction rate (this is further illustrated in Figs. 3]and|5}. The 
important point here is that the Monte Carlo sampling results 
in "low" and "high" rates for which the coverage probability 
can be quantified precisely. 

Although a low, high and median rate are useful quantities, 
they do not necessarily contain all the information on the rate 
probability density. For application of a reaction rate to nucle- 
osynthesis calculations, therefore, it is useful to approximate 
the rate probability density by a simple analytical approxima- 
tion. It was shown in 12211 that in most (but not all) cases 
the reaction rate probability density is well approximated by 
a lognormal distribution (equation|8}. The lognormal parame- 
ters /x and a can be found from the sampled total rates at each 
temperature according to 
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where _E[ln(y)] and V [ln(y)] denote the expectation value and 
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variance of the natural logarithm of the total rate, y, respec- 
tively. A useful measure of the applicability of a lognormal 
approximation to the actual sampled distribution is provided 
by the Anderson-Darling statistic^, which is calculated from 

" 2i — 1 

t AD = -n - 53 ( lnF (y*) + In [1 - F(y n+1 _i)] 

i=l U 

(13) 

where n is the number of samples, yi are the sampled reaction 
rates at a given temperature (arranged in ascending order), and 
F is the cumulative distribution of a standard normal func- 
tion (i.e., a Gaussian centred at zero). An A-D value greater 
than unity indicates a deviation from a lognormal distribution. 
However, it was found by Longland et al. 1I22I1 that the rate 
distribution does not visibly deviate from lognormal until A-D 
exceeds tAD ~ 30. The A-D statistic is presented in Tabs. I VII 
and IVHI along with the reaction rates at each temperature in 
order to provide a reference to the reader. 

C. Extrapolation of Experimental Reaction Rates to Higher 
Temperatures 

Experimental rates usually need to be extrapolated to high 
temperatures with the aid of theoretical models because res- 
onances are only measured up to some finite energy, E^ax- 
If the effective stellar burning energy window 113 811 extends 
above this energy, the rate calculated using the procedure out- 
lined above will become inaccurate. Statistical nuclear reac- 
tion models must, therefore, be used to extrapolate the exper- 
imental rates beyond this temperature. The method used here 
is described in detail in Ref. Dal . It uses the following strat- 
egy: (i) an effective thermonuclear energy range (ETER) is 
defined using the 8 th , 50 th , and 92 nd percentiles of the cumu- 
lative distribution of fractional reaction rates (i.e., the relative 
contribution of single resonances at temperature T divided by 
the total reaction rate at T); (ii) the temperature, T mata h, be- 
yond which the total rate must be extrapolated is estimated 
from 

£(T match ) + A£(T match ) = (14) 

where A£'(T matc h) is the width of the ETER calculated from 
the 8 th and 92 nd rate percentiles. We adopt the Hauser- 
Feshbach rates of Ref. 03911 for temperatures beyond T^tch, 
normalised to the experimental rate at T matc h. 

in. THE 22 NE+a REACTIONS 

A. General Aspects 

The 22 Ne(a,n) 25 Mg (Q an = -478.296(89) keV) 
and 22 Ne(a,7) 26 Mg (Q Q7 = 10614.787(33) keV) reac- 



3 The Anderson-Darling statistic t37ll is more useful than a x 2 statistic be- 
cause it does not require binning of the data. The latter usually results in a 
loss of information. 



tions are both important in s-process neutron production. 
While the 22 Ne(a,n) 25 Mg reaction produces neutrons, the 
22 Ne(a,7) 26 Mg reaction also influences the neutron flux by 
directly competing for available a-particles. The rates of 
both reactions will therefore be presented here. The centre - 
of-mass energy region of interest to the s-process amounts to 
E cm = 600 ± 300 keV, corresponding to excitation energies 
of E x = 10900 - 11500 keV in the 26 Mg compound nucleus. 
Note that only states of "natural" parity (i.e., + , l - , 2 + , etc.) 
can be populated via 22 Ne+a (because both target and pro- 
jectile have spin-parities of + ). 

Since the early 1980's, several direct measurements were 
performed of both reactions close to the energy region of 
interest 11191 140l 444j| . All of these measurements, with the 
exception of Ref. Il44ll . were made using gas targets at the 
Institut fur Strahlenphysik in Stuttgart, Germany [e.g., l45ll . 
The lowest energy resonance measured in those works is lo- 
cated at £y ab R3 830 keV, near the high energy end of the 
astrophysically important region. The structure of the 26 Mg 
compound nucleus near the a-particle and neutron thresholds 
has been investigated previously via neutron capture IkKh|47[| . 
scattering [ 481 - 15011 . photoexcitation Il5ll-l53ll . transfer 11441,154 - 
H^l, and photoneutron measurements ll57ll . In particular, the 
latter study observed the strong population of a 26 Mg level 
near E x = 11150 keV, with presumed quantum numbers of 

= 1~, corresponding to an expected low-energy resonance 
at E cm = 450 keV. It was believed to have been observed by 
Drotleff et al. JH and Harms et al. JH at E™> = 630 keV, 
but the presumed signal was later shown to be caused by back- 
ground from the 11 B(a,n) 14 N reaction. Nevertheless, the an- 
ticipated contribution from this low-energy resonance has sen- 
sitively influenced all past estimates of 22 Ne+a reaction rates. 
For example, it was shown by The et al. 15811 that it has a 
strong impact on s-process nucleosynthesis in massive stars. 
However, recent 26 Mg(7, 7) 26 Mg studies by Longland et al. 
l59ll demonstrated unambiguously that this particular level has 
unnatural parity (J 71 = 1 + ) and, therefore, cannot be popu- 
lated via a-particle capture on 22 Ne. 

Studies that provide new experimental information relevant 
to 22 Ne+a, obtained after the NACRE compilation was pub- 
lished [18], are summarised in table |T] The goal of the fol- 
lowing discussion is to consider all the available experimen- 
tal information for states in 26 Mg of interest to s-process nu- 
cleosynthesis and to assign these levels to corresponding res- 
onances in both 22 Ne(a,7) 26 Mg and 22 Ne(a,n) 25 Mg. This 
allows for an estimation of the partial and total resonance 
widths, resulting in more accurate 22 Ne+a reaction rates. 
A number of levels in 26 Mg near the a-particle and neu- 
tron thresholds have unknown spin-parities and partial widths. 
These levels have been disregarded in all previous reaction 
rate estimates. Since it is not known at present if any of these 
are natural parity states and, therefore, may be populated in 
22 Ne+a, they cannot be easily included in a Monte Carlo re- 
action rate analysis at present. Thus, our strategy is as follows: 
we will first derive 22 Ne+a Monte Carlo rates by excluding 
these levels of unknown spin-parities. Subsequently, we will 
investigate their impact on the total reaction rates under the 
extreme assumption that all of these levels possess natural 
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parity. As will be seen below, future measurements of these 
states are highly desirable. Throughout the following discus- 
sion, energies are presented in the centre of mass frame unless 
otherwise stated. 



B. Resonance Strengths 

Directly measured resonance strengths in the 
22 Ne(a,7) 26 Mg reaction in the energy range of 
E' r ab = 830 - 2040 keV are adopted from Ref. 
Direct measurements of resonances in the 22 Ne(a,n) 25 Mg 
reaction at ener gies of = 830 — 2040 keV are re- 
ported in Refs. HH HqI El |H Q. Note, however, that 
the 22 Ne(a,n) 25 Mg resonance strengths from the different 
measurements disagree by up to a factor of 5 (i.e., a deviation 
well outside the quoted uncertainties). Clearly, adopting a 
simple weighted average value would not account for the 
unknown systematic bias present in the data. To alleviate 
this problem, we adopt the method of Ref. [61], previously 
applied to account for unknown systematic uncertainties 
in neutron-lifetime measurements. This method follows a 
similar procedure for characterising unknown systematic 
uncertainties as that presented in Ref. 06211 . It assumes that 
all the reported strength values of a given resonance have the 
same, unknown, systematic error, a u , which can be summed 
in quadrature with the reported uncertainties. Hence for each 
reported uncertainty of data set i, er^, an inflated uncertainty, 
er', is obtained via 



hand, the inflated weighted average value is 017 = 1.4(3) x 
10~ 4 eV. We applied the inflated weighted average method to 
all resonances in the energy region E* = 830 — 1495 keV. 
Above this energy range, we used the (standard) weighted av- 
erage because the different data sets are in considerably better 
agreement. 

From the measured 22 Ne(a,n) 25 Mg and 22 Ne(a,7) 26 Mg 
strengths of a given resonance, the neutron and a-particle par- 
tial widths can be found if the 7-ray partial width can be esti- 
mated. This information allows for integrating the resonance 
cross section numerically, according to equation ((3), which is 
more reliable than adopting the narrow resonance approxima- 
tion, equation |0. Because of Coulomb barrier penetrability 
arguments, the neutron width is expected to dominate the total 
width of the resonances important for s-process nucleosyn- 
thesis (i.e., r„ sa T). Thus, in most (but not all) cases, the 
neutron width exceeds the a-particle width for a given state 
substantially and we can use the following approximations 
to determine the a-particle width from measured resonance 
strengths. For the 22 Ne(a,7) 26 Mg reaction, the a-particle 
partial width can be found by assuming a reasonable average 
value for the 7 -ray partial width of T 7 sa 3 eV [60]. We inves- 
tigated the effect of this choice on the reaction rates and the 
exact average value of T 7 was found to be relatively unim- 
portant. The a-particle partial width can then be found (for 
T n ps T) from 



W7 a7 = Ul 



Y Y 



Y a 7 T 7 7 Y r 



^7^7 r 

uj ieV 



(18) 



(15) 



For the 22 Ne(a,n) 25 Mg reaction, the a-particle partial width 
can be calculated from: 



From the inflated uncertainties, the weighted average of the 
resonance strengths, ojji, is obtained in the usual manner, 



UJJ ar 



(19) 



'2 



1 



(16) 



C. Spectroscopic Factors 



The unknown value of <j u is adjusted numerically until the 
reduced chi-squared, x 2 /v, becomes equal to unity. 



v 



1 



— Y 



(w7i - ojjY 



(17) 



where v is the degree of freedom (i.e., v = n—1, with n equal 
to the number of measurements). 

Application of this method has two consequences compared 
to calculating the weighted average of the reported resonance 
strength values: (i) the uncertainty of the resonance strength, 
<F, will be larger, reflecting the fact that the systematic shift 
in the data is of unknown nature; and (ii) strength values 
with small reported uncertainties will carry less weight. Con- 
sider as an example the lowest-lying observed resonance in 
the 22 Ne(a,n) 25 Mg reaction, located at = 831 keV. The 
measurements reported in Refs. lfl9l WH R3I R4I1 yield for the 
resonance strength a (standard) weighted average of uij = 
1.2(1) x 10~ 4 eV, with \ 2 jv = 2.9, indicating poor agree- 
ment between the individual measurements. On the other 



Alpha-particle spectroscopic factors for levels near the a- 
particle and neutron thresholds in 26 Mg have been obtained 
from 22 Ne( 6 Li,d) 26 Mg transfer studies by Refs. HH and 0]. 
The spectroscopic factors derived from the ( 6 Li,d) transfer 
data are important because they allow for an estimate of the 
a-particle partial width, T a , of 22 Ne7a resonances via equa- 
tions ©, ©, and ©. 

Numerous studies have shown that a-transfer measure- 
ments are very useful for measuring relative spectroscopic 
factors, but are not sufficiently accurate for predicting abso- 
lute values. For this reason, the measured spectroscopic fac- 
tors are frequently scaled relative to resonances with well- 
known partial widths (note that this is an approximation 
equivalent to assuming that 9 2 = S in Sec. [ilj. For example, 
Giesen et al. 114411 scaled their spectroscopic factors relative to 
the E l * h = 831 keV ( J n = 2+ ) resonance in 22 Ne(a,n) 25 Mg. 
Our best value for the (a,n) resonance strength is LJj an = 
1.4(3) x 10~ 4 eV (see Tab. lTVT i. Since for this low-energy res- 
onance it can be safely assumed that T pa T n , a spectroscopic 
factor of Sa' n ^ = 0.98 is obtained from equations © and (0. 
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Reference 


Reaction Studied 


Comments 


TflpQPr ft nl r 1 Ql 

Koehler [60J 
Ugalde et al. [56] 
Longland et al. [59] 
DeBoer et al. [2i] 


22 Ne<7> n^ b MsJ 

1 i L^l f_X ■ 11 / i-V±il 

na, Mg(n, 7 ) 
22 Ne( 6 Li,d) 26 Mg 
26 Mg(7, 7 ') 26 Mg 
26 Mg(7, 7 ') 26 Mg 


Resonances hetween E VjS> — 570 keV and E ub — 1450 keV 

E x ,J n ,T y , T„ for states corresponding to E^ h = 570 keV to E^ h = 1000 keV 

J w and 5 Q for two states below neutron threshold 

E x and J" for four resonances corresponding to E ] * b = 38 keV to E^ h = 636 keV 
T 7 for four resonances corresponding to E^° = 38 keV to E^ h — 636 keV 



TABLE I. New information relevant to the 22 Ne+« reaction rates that has become available since the NACRE compilation was published 1 18]. 



Surprisingly, this value is a factor of 27 larger than the transfer 

value extracted by Ref. H, S"i 6Ll ' d) = 0.037. Renormalisa- 
tion of all measured ( 6 Li,d) spectroscopic factors to the (a,n) 
spectroscopic factor of the £y ab = 831 keV resonance results 
in the values shown in green in Fig. [p It is certainly remark- 
able that all levels observed by Ref. 14411 should have dimen- 
sionless reduced a-particle widths far larger in value than the 
Porter- Thomas prediction. Additionally, several of these lev- 
els exhibit dimensionless reduced widths near or exceeding 
the Wigner limit, even if one accounts for the difference be- 
tween S a and 6 2 a (Sec. Ill Al ). 

However, there is no compelling reason why the £y ab = 
831 keV resonance should be singled out for the normalisation 
procedure, other than it being the lowest-lying observed res- 
onance. For example, one may consider another well-known 
resonance, located at £y ab = 1434 keV (J 7 * = 2+ ). From its 
measured 22 Ne(a,n) 25 Mg resonance strength, an a-particle 
spectroscopic factor of Sa* ,n ^ = 0.27 is obtained. The a- 
particle transfer value for the corresponding levels, measured 

by Ref. i44ll . amounts to Sa L> ' d ^ = 0.11. These two values 
differ by a factor of 2.5, and thus are much closer in agreement 
than the results for the E^ h = 831 keV resonance. Normali- 
sation of all measured relative ( 6 Li,d) spectroscopic factors to 
the (a,n) spectroscopic factor of the E x f° = 1434 keV reso- 
nance results in the values shown in blue in Fig.Q] It is evident 
that these normalised values are in far better agreement with 
the Porter-Thomas distribution than the results obtained when 
scaling spectroscopic factors relative to the Ey ab = 831 keV 
resonance. In addition, by using the _Ej? b = 1434 keV res- 
onance normalisation, all of the resulting dimensionless re- 
duced widths now have values less than the Wigner limit, 
making them more believable. Since we feel it is more rea- 
sonable to scale the ( 6 Li,d) spectroscopic factors using the 
El ah = 1434 keV resonance instead of the E l r ab = 831 keV 
resonance, we adopt the reduced widths shown in blue in 
Fig.[T]for calculating the 22 Ne+a rates. Note that the relative 
spectroscopic factors obtained by Ref. II44|| have been used at 
face value by Refs. lfl2l l60ll in their reaction rate calculations. 
Clearly, this issue needs to be resolved in future work. 



IV. INFORMATION ON SPECIFIC 26 MG LEVELS 

E x = 10693 keV (E' r ab = 92 keV; J 71 = 4+ ). An ex- 
cited state near this energy has been observed by Glatz et al. 
JH at E a =10695f2) keV, Giesen et al. H at E a; =10694(20) 
keV, and Moss H at £^=10689(3) keV. A weighted aver- 



age of these excitation energies is used in the present work. 
The J 71 " = 4 + assignment was made by considering the de- 
cay scheme of this state as observed by Ref. [54] and that the 
state most likely has natural parity. The a-particle spectro- 
scopic factor for this state from Ref. [44], after normalisation, 
is S a = 0.059. 

E x = 10806 keV (E' r ab = 226 keV; J"" = 1" ). This 
state was seen previously in 22 Ne( 6 Li,d) 26 Mg measure- 
ments by Ugalde et al. OS at £^=10808 (20) keV, and in 
25 Mg(n t ,7) 26 Mg measurements (thermal neutron capture) by 
Walkiewicz et al. JH at E x =10805.9 (4) keV. A recent exper- 
iment assigned a spin-parity of J T = 1~ ll59ll . The adopted 
excitation energy is the weighted average of these results. The 
a-particle spectroscopic factor from Ref. 115 611 . after normali- 
sation, amounts to S a = 0.048. 

E x = 10943 keV (E' r ab = 388 keV; 3 n = (5" - 7 - )). 
An excited state at this energy has been observed by Glatz 
et al. l54Tl . The observed decay scheme restricts the quan- 
tum numbers, using the dipole-or-E2 rule of Ref. 0641, to 
J n = (5 ± - 7~). A 22 Ne( 6 Li,d) 26 Mg transfer measurement 
populated a state at £^=10953 (25) keV, but did not obtain the 
quantum numbers other than to report that it most likely had 
natural parity ll56ll . The combined quantum number assign- 
ment is therefore J 71 = (5~ — 7~). This state is treated here as 
part of a doublet with the E x = 10949 keV state. Note that in 
the reaction rate calculations of Karakas et al. IU2I1 . this state 
was incorrectly assigned spin-parity values of J x = 2 + , 3~. 

E x = 10949 keV (E' r ab = 395 keV; J n = 1~ ). This 
state has been observed previously in 26 Mg(p,p') 26 Mg mea- 
surements by Moss dm at E x = 10950(3) keV. The (p,p') 
measurements suggest a J" = 1~ assignment, which agrees 
with the 26 Mg(7, 7) 26 Mg result of Longland et al. It 
is unclear whether Ugalde et al. B56I1 observed this state or 
the one at Ej;=10943 keV. Therefore, in the present analy- 
sis, the normalised spectroscopic factor of S = 7 x 10~ 3 re- 
ported in Ref. ll56ll is treated as an upper limit for both states 
at E x = 10943 keV and E x = 10949 keV. 

E x = 11112 keV (E lab = 587 keV; J T = 2+ ). The state 
is located above the neutron threshold. It has been observed 
previously in a 25 Mg(n,7) 26 Mg experiment [47, 60] and was 
assigned a spin-parity of J 77 = 2+ . 

E x = 11154 keV (E lab = 637 keV; J ff = 1+ ). A state 
at this energy has been observed by Fagg [48], Weigmann 
et al. |H, Crawley et al. HI, Yasue et al. J5J], Koehler ll60ll . 
Tamii et al. 15011 . and Schwengner et al. l53ll . Additionally, 
an excited state at this energy was strongly populated by the 
photoneutron experiment of Berman et al. 15711 . who predicted 
a J" = 1~ assignment. As a result of this prediction, sev- 
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FIG. 1. (Colour online) Dimensionless reduced a-particle widths of unbound states from Ref. [36], and references therein [see also. l22h . Also 
plotted is the Porter-Thomas distribution that best fits these data at small values. It is apparent from the figure that states with large a-particle 
spectroscopic factors are not represented by the Porter-Thomas distribution. These levels most likely have an a-particle cluster structure and 
would be populated preferentially in transfer measurements, such as the ( 6 Li,d) measurements of Refs. 1 44] and l56ll . Also shown in green and 
blue are the normalised spectroscopic factors measured by Ref. 1441 . The values normalised using the E^ h = 1434 keV resonance are shown 
in blue, while those normalised to the _Ej? b = 831 keV resonance are shown in green. Clearly, the normalisations are vastly different, and the 
spectroscopic factors obtained using the Er h = 831 keV resonance as a normalisation reference appear to be too high, as shown in the figure 
inset, which displays the same information but on an expanded scale. See text for more detail. 



eral studies have searched for a resonance corresponding to 
this energy |fl9l I4ll - l44l l56ll . Of these studies, a presumed res- 



onance was reported by Refs. ll4l"ll42 
caused by beam induced background 



Mg(7, 7) Mg experiment (5 



, but later proven to be 
4311 . Recently, however, 
showed unambiguously 



that the spin-parity of this state amounts to J 11 = 1 + and, 
therefore, cannot contribute to the 22 Ne+a reactions rates. 
A more detailed discussion of this state is presented in sec- 
tioning] 

E x =11163-11326 keV (E' r ab = 648 840 keV). Excita- 
tion energies were taken as weighted averages of Moss ll49ll . 
Glatz et al. J3, and Koehler [60]. Quantum numbers, neu- 
tron and 7-ray partial widths were all adopted from Ref. ll60ll . 
Since no a-particle partial widths have been measured for 
these states, upper limits were derived either from the data 
presented by Ref. 14011 . or adopted from the maximum theo- 
retically allowed values, depending on which was smaller. 

E x = 11318 keV (E' r ab = 831 keV; J" = 2+ ). Koehler 
ll60ll argued that this state cannot correspond to both the res- 
onance observed by Jaeger et al. B19I1 at E r =832 (2) keV in 
22 Ne(a,n) 25 Mg and by Wolke et al. 10] at E l r ab = 828(5) 



keV in 22 Ne(a,7) 26 Mg, because the implied value of T 7 = 
76 eV would be far larger than the average 7-ray partial 
width (T 7 = 3 eV) in this energy range. However, this 
conclusion is questionable considering the large uncertainty, 
T 7 = 76(53) eV, when the 7-ray partial width is derived from 
the measured values of wy ai , ojj an , and T. Clearly, the de- 
viation from the average in this energy range amounts to only 

1.4(7. 

Since it cannot be decided at present if the (a,n) and (a, 7) 
resonances correspond to the same 26 Mg level or not, the par- 
tial widths cannot be derived unambiguously from the mea- 
sured resonance strengths and total width. Thus we assumed 
that the (a,n) and (a, 7) resonances are "narrow", i.e., we em- 
ployed equation (O instead of equation (|3]l in our rate cal- 
culations. The strength reported by Ref. fl40n is used for the 
22 Ne(a,7) 26 Mg resonance, while the inflated weighted aver- 
age (see section lHIBI ) is adopted for the 22 Ne(a,n) 25 Mg reso- 
nance, resulting in a strength of wyr a n -\ = 1-4(3) x 10~ 4 eV. 

E x =11328-11425 keV (E' r ab = 843 - 957 keV). Excita- 
tion energies, quantum numbers, neutron, and 7-ray widths 
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for these levels are adopted from Refs. ll47ll . and ll60ll . No a- 
particle widths have been measured for these states, and thus 
upper limits have been adopted from either the data presented 
by Refs. B40I1 and 11911 . or from the maximum theoretically 
allowed values, depending on which was smaller. 

E x > 11441 keV (E|? b > 976 keV). Resonances corre- 
sponding to excited states above E x = 11441 keV have been 
measured directly [ 19, 40-44]. In order to take the widths of 
wide resonances into account, the neutron and 7-ray partial 
widths (and quantum numbers) measured by Refs. Il47ll and 
ll60ll have been used when available. The inflated weighted 
average method (see section lHIBI i is used to combine the dif- 
ferent 22 Ne(a,n) 25 Mg strengths for resonances below E^ h — 
1434 keV, while standard weighted averages are used above 
this energy. Since the 22 Ne(a,7) 26 Mg resonances measured 
in Ref. [40] cannot be assigned unambiguously to correspond- 
ing 22 Ne(a,n) 25 Mg resonances, all of these resonances were 
treated as independent and narrow. The quantum numbers of 
22 Ne(a,n) 25 Mg resonances located above E x * h = 1530 keV 
are adopted from Ref. ll40ll when not available otherwise. 



V. REACTION RATES FOR 22 NE+q 

The resonance properties used to calculate the rates for both 
the 22 Ne(a,7) 26 Mg and the 22 Ne(a,n) 25 Mg reactions are pre- 
sented in Tabs. [II]— [V] For more detailed information on level 
properties, see Ref. 113011 - Separate tables are used to list reso- 
nances with measured partial widths and those which possess 
only an upper limit for the a-particle width but have known 
neutron and 7-ray widths. 

The matching temperature, T matc h, (see Sec. Ill CI ) for 
both the 22 Ne(o!,7) 26 Mg and 22 Ne(a,n) 25 Mg reactions, be- 
yond which the rates are estimated by normalising Hauser- 
Feshbach predictions to experimental rates, amounts to T = 
1.33 GK i.e., well above the temperatures relevant for the s- 
process during He-burning (T = 0.01 — 0.3 GK). 

Monte Carlo reaction rates for the 22 Ne(a,7) 26 Mg and 
22 Ne(a,n) 25 Mg reactions are presented in Tabs. [VI] and 
IVIII respectively. The median, low, and high rates are 
shown alongside the lognormal parameters and the Anderson- 
Darling statistic described in Sec. Ill Bl The Monte Carlo reac- 
tion rate probability density functions are displayed in Figs. [2] 
and[3]as red histograms. The solid black lines indicate the log- 
normal approximation, calculated with the lognormal param- 
eters, /1 and a, listed in columns 5, 6, 10, and 11 of Tabs. [VI] 
andED 

In order to emphasise that our low and high rates, obtained 
for a coverage probability of 68% (see section HTBT i. do not 
represent sharp boundaries, we show the (0,7) and (a,n) re- 
action rates, normalised to the respective recommended (me- 
dian) values, as colour contours in Figs. |4]and[5] The thick 
and thin solid lines represent coverage probabilities of 68% 
and 95%, respectively. The three dashed lines show the pre- 
viously reported rates (Angulo et al. iQjl for 22 Ne(a,7) 26 Mg 
and Jaeger et al. Il9ll for 22 Ne(a,n) 25 Mg), normalised to our 
recommended rate. Our calculations of the relative resonance 
contributions to the total (a, 7) and (a,n) reaction rates show 



that, at temperatures most relevant to the s-process, reso- 
nances including and below the E x f° = 831 keV resonance 
are the most important. Future experimental efforts should, 
therefore, be concentrated on studying resonances in the exci- 
tation energy region near the neutron threshold. 

For the 22 Ne(a,7) 26 Mg reaction, the present rates deviate 
significantly from the results of Ref. [ 18], by factors of 2-100. 
The differences are caused by: (i) a different treatment of par- 
tial widths; in Ref. lfl8ll the rates were found from numerical 
integration by assuming upper limit values (T = 4 — 10 keV) 
for the total widths, whereas in the present work total widths 
have been adopted from measured values; (ii) our improved 
treatment of upper limits for reduced a-particle widths (i.e., 
sampling over a Porter- Thomas distribution; see section ITTBl ; 
and (iii) the fact that new nuclear data became available since 
1999 (see Tab. [I]). The combined effect of these improvements 
results in a factor of 5 reduction in reaction rate uncertainties 
in the He-buming temperature region. 

As already noted in section IHI Al a number of excited 
states near the a-particle and neutron thresholds in 26 Mg have 
been observed by additional inelastic proton scattering exper- 
iments B491 pOl I52II . However, their spins and parities have 
not been determined. In particular, it is not known at present 
if these levels possess natural parity and thereby may be pop- 
ulated in the 22 Ne+a reactions. In order to investigate the 
maximum impact of these states with unknown J 71 ' values on 
the 22 Ne(a,7) 26 Mg reaction rate, we performed a test by as- 
suming that all of these levels possess natural parity and by 
adopting upper limit a-particle spectroscopic factors from that 
data of Refs. 114-411 and 15611 . The results show that these states 
can increase the 22 Ne(a,7) 26 Mg reaction rate by up to a factor 
of 30 at temperatures between T9 = 0.1 and 0.2 GK. 

For the 22 Ne(a,n) 25 Mg reaction there is better agreement 
between previous and new rates. The present rates are slightly 
higher (up to a factor of 2) than those calculated by Ref. 111911 . 
The two main reasons for the difference are: (i) we used in- 
flated weighted averages of the reported resonance strengths 
from different measurements (see section IIHBI ; and (ii) ex- 
cluded the contribution of a presumed E^ h = 630 keV reso- 
nance, because the level at E x — 11154 keV has been shown 
to possess unnatural parity 05911 . 

We would like to emphasise that the observed (a,n) and 
(a,7) resonances near E^ h = 830 keV introduce another sys- 
tematic uncertainty that we have not accounted for. Recall 
that we treated these two resonances as independent and nar- 
row (sectionlTVli. On the other hand, if they correspond to the 
same level in 26 Mg, the partial widths could be derived from 
the measured resonance strengths. In that case, the resonance 
turns out to be relatively broad, resulting in a significant con- 
tribution of the resonance tail to the total reaction rate. Tests 
show that the resulting reaction rates near T w 0.3 GK could 
increase by roughly a factor of 5. 

The ratio of 22 Ne(a,n) 25 Mg to 22 Ne(a,7) 26 Mg reaction 
rates is shown in figure [6] Note that these rates are not in- 
dependent since, for example, the same values of a-particle 
partial widths enter in both rate calculations if an (a, 7) and 
(a, n) resonance corresponds to the same 26 Mg level. Thus 
the uncertainties shown in Fig.[6]are somewhat overestimated. 
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T9 = 0.06 
A-D= 124.5 
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Reaction rate (cm 3 mol 1 s 1 ) 

FIG. 2. (Colour online) Reaction rate probability densities for the 22 Ne(«,7) 26 Mg reaction at various stellar temperatures. In each panel, the 
red histogram represents the Monte Carlo results, while the solid line shows the lognormal approximation. Note that the solid line is not a fit 
to the histogram, but was calculated from the lognormal parameters fi and a (table |V0 , which in turn were determined from equation (I12t . It 
is apparent that the lognormal approximation to the reaction rates holds in the temperature range of the s-process (near 0.3 GK). 
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FIG. 3. (Colour online) Reaction rate probability densities for the 22 Ne(a,n) 25 Mg reaction. See caption to Fig.f2] 
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a Average value from Ref. 16011 

b From 22 Ne(o,n) 25 Mg measurements (see equation H9\ ) 

c Detailed discussion on quantum number assignments can be found in section lTtlBI 

TABLE II. Resonances of 22 Ne(a,7) 26 Mg with known a-particle partial widths or resonance strengths. Total widths are from Ref. 14011 for 
resonances above E^ h = 1533 keV. For lower-lying resonances, total widths are adopted from Ref. Ill9ll and Ref. 16011 . Ambiguous spin- 
parities (i.e., those not based on strong arguments) are placed in parentheses, according to the guidelines in Ref [64]. The last column, labelled 
"Int" indicates those resonances for which sufficient information is available in order to integrate their reaction rate contribution numerically, 
according to equation l[3}. 



Nevertheless, it is instructive to compare the present ratios, 
shown in black, to those from previous work 111 8l Il9ll . dis- 
played in red. It can be seen that the present ratio is sig- 
nificantly larger than previous results and, consequently, we 
predict that more neutrons will be produced per captured a- 
particle. 



VI. ASTROPHYSICAL IMPLICATIONS 



rates will be correlated, it is difficult to account for these cor- 
relations since the 22 Ne(a,7) 26 Mg rate includes resonances 
below the neutron threshold and since some 26 Mg levels con- 
tribute more to one reaction channel than the other. For these 
reasons, we have chosen in the present study to explore con- 
servatively the impact of the largest reaction rate variations. 
These nucleosynthesis calculations are performed separately 
for massive stars and AGB stars. 



A. Models 



1. Massive Star Models 



In order to explore how the current 22 Ne+a reaction 
rates affect s-process nucleosynthesis, two kind of calcula- 
tions are presented here. The first compares final abundance 
yields from post-processing models upon changing the rec- 
ommended 22 Ne+a reaction rates from previously published 
results to those presented in this paper. The second calcula- 
tion estimates the variations in s-process nucleosynthesis aris- 
ing from uncertainties in the present 22 Ne+a reaction rates. 
These can then be compared with abundance variations aris- 
ing from the literature rates. 

In order to take the uncertainties into account, three sets 
of calculations were performed: (i) recommended rates for 
both 22 Ne(a,n) 25 Mg and 22 Ne(a,7) 26 Mg reactions, (ii) low 
22 Ne(a,n) 25 Mg rate and high 22 Ne(a,7) 26 Mg rate, and (iii) 
high 22 Ne(a,n) 25 Mg rate and low 22 Ne(a,7) 26 Mg rate. Al- 
though in reality the 22 Ne(a,7) 26 Mg and 22 Ne(a,n) 25 Mg 



A single zone temperature-density profile has been used 
to study the effects of the 22 Ne+a reaction rates on nu- 
cleosynthesis during the core helium burning stage in mas- 
sive stars. The temperature-density profile and initial abun- 
dances used in the present study are for a 25M star and 
have been used previously in Refs. lf24l l58ll . The most abun- 
dant isotopes at the onset of helium burning are (in mass frac- 
tions, X): 4 He (X a = 0.982), 14 N (Xi4 N = 0.0 1 22), 20 Ne 
(X 2 o Ne = 0.0016), and 60 Fe (X 6 o Fe = 0.00117). During most 
of the core helium burning phase, the temperature and density 
(T rj 100 - 250 MK and p m 1000 - 2000 g/cm 3 , respec- 
tively) are not high enough for the 22 Ne(a,n) 25 Mg neutron 
source to produce a significant number of neutrons. How- 
ever, towards the end of this phase the temperatures become 
high enough for efficient neutron production. The exact time 
at which the 22 Ne(a,n) 25 Mg reaction starts to occur not only 
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TABLE III. Properties of Unobserved Resonances in Ne(a,7) Mg . For these resonances, only upper limits of the resonance strength and/or 
the a-particle spectroscopic factor are available at present. The 7-ray and neutron partial widths are taken from the R-matrix fit of Ref. |60], 
Quantum numbers for states below E x — 11163 keV are discussed in Sec. lIVI All other quantum numbers are adopted from Ref. [ 60] . When 
a range of quantum numbers is allowed, the upper limit of the a-particle width is calculated assuming the lowest possible orbital angular 
momentum transfer. The upper limit a-particle spectroscopic factors adopted (Sc.vl) are also listed for completeness. 



affects the number of neutrons produced during core helium 
burning, but also the amount of 22 Ne remaining that can be 
processed later during the carbon shell burning phase. Al- 
though not studied here, the s-process is also expected to be 
active during shell carbon burning. 

The nucleosynthesis study was performed with a 583 nu- 
cleus s-process network that extends up to molybdenum. Re- 
action rates (other than the 22 Ne+a rates) were adopted from 
the star lib library ||65ll . The star lib library incorpo- 
rates a compilation of recently evaluated experimental Monte 
Carlo reaction rates in tabular format on a grid of 60 temper- 
atures from 1 MK to 10 GK. Tabulated are the temperature, 
the reaction rate, and the factor uncertainty, which is closely 
related to the lognormal parameter, er, in Ref. l22tl . 

2. AGB Star Models 

The AGB nucleosynthesis tests are performed on a 5.5 M Q , 
Z = 0.0001 model star, detailed in Ref. Q. This model 
was chosen because it experiences many thermal pulses (77 in 
total) during the AGB phase, where 69 of those He-shell in- 
stabilities reach peak temperatures of 0.30 GK or higher (with 
temperatures of 0.35 GK for 50 thermal pulses). One com- 
plication arises from disentangling the effects of the 22 Ne+a 
rates and those of proton-capture nucleosynthesis at the base 
of the convective envelope (hot bottom burning, HBB). In our 



model, the base of the envelope reaches peak temperatures 
of 98 MK, easily hot enough for activation of the NeNe and 
MgAl proton-burning chains. The main results were reduc- 
tions in the envelope 24 Mg and 25 Mg abundances, and in- 
creases in 26 Mg, 26 Al, and 27 A1. This means that the He- 
intershell preceding a pulse contains a non-solar Mg isotopic 
composition that is enriched in 26 Mg. 

The post-processing nucleosynthesis used for the AGB star 
models has previously been described in detail by Ref. [e.g., 
16611 . This code needs as input from the stellar evolution code 
variables such as temperature, density, and convective bound- 
aries as a function of time and mass fraction. The code then 
traces the abundance changes as a function of mass and time 
using a nuclear network containing 172 species (from neu- 
trons to sulphur, and then from iron to molybdenum) and 
assuming time-dependent diffusive mixing for all convective 
zones 1 67]. Although this network does not contain species 
of the main s-process above A 1=3 100, their production is es- 
timated by the inclusion of an extra isotope (the "g particle", 
counting neutron captures beyond our network). The reac- 
tion rates used in the nuclear network are mostly taken from 
the JINA reaclib database [68], with the exception of the 
22 Ne+a rates adopted from the present work. Some modifica- 
tions were made to the JINA reaclib library including the 
removal of the 96 Zr decay rate (since this is an essentially sta- 
ble isotope with a half-life of t 1 / 2 > 10 19 years), and the in- 
clusion of the ground and isomeric states in 85 Kr. This is done 
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il A detailed discussion of quantum number assignment can be found in the text. 

b Calculated using equation (19). 

c Average value from Ref. I60ll . 

d Assuming V is dominated by T n (see section ITlI At . 

TABLE IV. Resonances in 22 Ne(a,n) 25 Mg with known a-particle partial widths or resonance strengths. When a range of quantum numbers is 
present, the one used for calculating the reaction rates is presented in bold. 



because 50% of the neutron flux from n+ 8 Kr proceeds to the 
ground state of 85 Kr (t 1/2 = 3934.4 days) and the other 50% 
goes to the isomeric state (t = 4.480 hours). The inclusion of 
both 85 Kr states is essential for Rb abundance predictions in 
AGB nucleosynthesis models [see discussion inM l69ll . 



B. Results 



The effects of our new rates on the nucleosynthesis in com- 
parison to using the results obtained in the literature are shown 
in Fig. [7] The improvements in abundance predictions for 
the two stellar environments are shown in Fig. [8] The most 
up-to-date previously published rates for the 22 Ne(ct,7) 26 Mg 
and 22 Ne(a,n) 25 Mg reactions are from Refs. IU8I1 and IU9I1 . 
respectively. The effects are markedly different for the two s- 
process environments, hence they will be discussed separately 
in the following. 



1. Massive Stars 



The recommended Ne(ct,n) J Mg reaction rate has not 
changed significantly in the present analysis. Consequently, 
we do not expect the final 25 Mg abundance to change. The fi- 
nal 26 Mg abundance, on the other hand, changes significantly 
by roughly a factor of three. The abundance changes in nu- 
clei heavier than iron are smaller, with the largest abundance 
increases occurring near 64 Ni. The increased destruction of 
isotopes already present in the star is also apparent for the p- 
nuclides 74 Se, 78 Kr, 84 Sr, and 93 Nb. These results indicate 
that with the reduced 22 Ne(a,7) 26 Mg rate, more neutrons are 
produced per 22 Ne+a reaction. Rather than extending the 
reach of the weak s-process component (i.e., synthesis of more 
massive nuclei), this flux increase affects branchings in the s- 
process path close to the iron peak. A wider range of inter- 
mediate mass nuclei are therefore produced. Fig. Qalso illus- 
trates that the 22 Ne+a rates not only affect the abundances 
of traditional s-process nuclides, but also the abundances of 
nuclei below the iron peak that act as poisons. An example is 
25 Mg, which produces 26 Mg through the 25 Mg(n,7) 26 Mg re- 
action. With a higher flux of available neutrons, this neutron 
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TABLE V. Properties of Unobserved Resonances in 22 Ne(a,n) 25 Mg . For these resonances, only upper limits of the resonance strength and/or 
the a-particle spectroscopic factor can be derived. Quantum numbers, 7-ray and neutron partial widths are taken from the R-matrix fit of Ref. 
l6Ch . Resonance energies represent a weighted average of values adopted from Refs. I4^7ll49ll54ll60ll . 



poison reaction occurs more frequently, effectively lessening 
the impact of the increased neutron flux on s-process nucle- 
osynthesis. 

Uncertainties in s-process nucleosynthesis in massive stars 
arising from uncertainties in the 22 Ne+a reaction rates are 
shown in Fig. [5] where the thin (red) bars show uncertainties 
arising from the old rates, and thicker (black) bars show those 
from the new rates. In particular, large reductions are notice- 
able for 26 Mg, where the current yield uncertainty amounts to 
around 50% in contrast to the previous factor of 5. Uncertain- 
ties in weak s-process nucleosynthesis have also undergone 
significant improvements, especially for species that can only 
be destroyed, but not created, by neutron captures. An ex- 
ample of this is the nucleus 58 Ni whose yield uncertainty has 
been reduced from a factor of five to just 50%. It is important 
to note here that, although the Monte-Carlo reaction rates do 
take into account systematic uncertainties, it is difficult to ac- 
count for ambiguities in the data, for example, the open ques- 
tion of whether or not the E~ b = 830 keV resonance is a 
doublet. Clearly, more measurements are needed. 



2. AGB Stars 

Nucleosynthesis yields from our low metallicity AGB star 
models show a very different pattern to those of the massive 
star study. For AGB stars, the effect on lighter elements is 
reduced in comparison to massive stars, with higher mass s- 
process elements revealing the largest changes. This weight- 
ing toward higher mass s-process elements is caused by our 
choice of using a low metallicity model. At low metallicity, 



the neutron/Fe seed ratio is much higher meaning that there 
is a higher production of higher atomic mass nuclei (e.g., see 
discussion in Refs. |H|,|7(3]). Nuclei towards the upper end of 
our network are produced up to a factor of 2 more than before, 
with the 'g' particle representing nuclei beyond our network 
capturing over 70% more neutrons. In low metallicity AGB 
stars, therefore, the 22 Ne+a reactions can be expected to pro- 
duce more high-mass s-process elements, while leaving the 
low-mass s-process below A ps 80 largely unaffected. 

Uncertainties in s-process nucleosynthesis have been, as in 
massive stars, dramatically improved with our new rates. The 
previous abundance uncertainties were approximately a fac- 
tor of 10, while the present uncertainties amount to less than 
a factor of 2. The present uncertainties in the rates affect the 
lower masses from Asa25 to Ass35 more than the s-process 
abundances. The ratio of 26 Mg and 25 Mg is still uncertain by 
approximately 20%, whereas it was previously around 80% 
(note that Ref. U2] found 26 Mg/ 25 Mg ratio uncertainties of 
60%). Rubidium and zirconium isotopes have undergone 
yield uncertainty improvements by a factor of about two. For 
the s-nuclide 96 Mo, the uncertainty has been reduced from a 
factor of 4 to a factor of 2 with our present results. 

The new 22 Ne+a reaction rates presented here should also 
be tested with low-mass AGB star models (M < 3M ). In 
lower mass AGB stars, while the 13 C(a,n) le O reaction is the 
main neutron source active between thermal pulses, activation 
of the 22 Ne+a reactions during a convective thermal pulse 
can have a significant effect on branchings in the s-process 
path. 
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TABLE VI. Monte Carlo reaction rates for the Ne(a,7) Mg reaction. Shown are the low, median, and high rates, corresponding to the 16th, 
50th, and 84th percentiles of the Monte Carlo probability density distributions. Also shown are the parameters (/i and <r) of the lognormal 
approximation to the actual Monte Carlo probability density. See Ref. l22ll for details. The rate values shown in parentheses indicate the 
temperatures (T > T matc h = 1.33 GK) for which Hauser-Feshbach rates, normalised to experimental results, are adopted (see section ITTCl . 
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TABLE VII. Monte Carlo reaction rates for the 22 Ne(a,n) 2j Mg reaction. Shown are the low, median, and high rates, corresponding to the 16th, 
50th, and 84th percentiles of the Monte Carlo probability density distributions. Also shown are the parameters (/i and <r) of the lognormal 
approximation to the actual Monte Carlo probability density. See Ref. (22h for details. The rate values shown in parentheses indicate the 
temperatures (T > T matc h = 1-33 GK) for which Hauser-Feshbach rates, normalised to experimental results, are adopted (see section ITTCl . 
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FIG. 4. (Colour online) The uncertainty bands for the 22 Ne(a,7) 26 Mg reaction. The uncertainties are the result of upper limit resonance 
contributions and of resonance strength uncertainties. The colour-densities represent the present reaction rate probability densities normalised 
to our recommended rate. The thick and thin black lines represent the 68% and 95% uncertainties, respectively. The dashed blue lines represent 
the literature rates from Ref. [18], with the thick and thin lines denoting the recommended rate and rate limits, respectively, normalised to our 
recommended rate. Values below unity (dotted line) indicate that the rates are lower than the present recommended rate. The relevant 
temperatures for helium- and carbon shell-burning have been added as red bars with the labels "He" and "C", respectively. 



VII. CONCLUSIONS 

Both the 22 Ne(a,n) 25 Mg and the 22 Ne(a,7) 26 Mg reactions 
influence the neutron flux available to the s-process in mas- 
sive stars and AGB stars. Uncertainties in the rates, there- 
fore, lead to large uncertainties in s-process nucleosynthesis. 
In this paper, we have estimated greatly improved 22 Ne+a 
reaction rates, based on newly available experimental infor- 
mation published since the works of Ref s. IU9ll and lfl8ll . and 



by applying a sophisticated rate computational method l22ll . 
Subsequently, we explored the astrophysical consequences for 
massive stars and for AGB stars. 

In massive stars, simple one zone models of core helium- 
burning were utilised to determine the influence of the new 
rates on the weak component of the s-process. The most im- 
portant result of our study is a significant reduction of nucle- 
osynthesis uncertainties. The yield uncertainty has been re- 
duced by between a factor of 5 and 10 across the s-process 
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FIG. 5. (Colour online) The probability densities for the 22 Ne(a,n) 25 Mg reaction in comparison to those presented by Ref. [19]. See caption 
of Fig.|4]for an explanation. 



mass region considered here (A < 100). For example, the 
yields of key isotopes, 26 Mg and 70 Zn, have uncertainty re- 
duction factors of about 5 and 10, respectively. When com- 
paring abundances obtained from our new recommended rates 
with those derived from previous recommended rates, the fi- 
nal yield of 26 Mg is found to have been reduced by roughly 
a factor of three, while s-process isotopes were affected only 



marginally. However, s-process nucleosynthesis is more con- 
centrated around the iron peak when using the new reaction 
rates. This relative insensitivity to changes in neutron flux is 
partially caused by captures on the neutron poisons 12 C, le O, 
and 25 Mg, which are present in large quantities. 

In our AGB star models, the final abundance uncertainties 
have also been improved significantly with the new rates, with 
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FIG. 6. (Colour online) Uncertainty bands of the reaction rate ratio, NA{o~v)( a , n )/NA{o~v)( ,,~ l ). The solid (black) lines represent the 
present reaction rate ratio, while the dashed (blue) lines represent the ratio of rates from Ref. |jl9ll for 22 Ne(a,n) 25 Mg, and Ref. fiUl . for 
22 Ne(a,7) 26 Mg. The recommended ratio (the center line in each set) was calculated by dividing the recommended 22 Ne(a,n) 25 Mg reac- 
tion rate by that of the 22 Ne(a,7) 26 Mg reaction at each temperature. To obtain the uncertainty bands for the rate ratio, the high rate for 
22 Ne(Q, 7 ) 26 Mg was divided by the low rate for 22 Ne(a,n) 25 Mg, and vice versa. Values greater than unity indicate that more neutrons (and 
25 Mg) are produced than 7-rays (and 26 Mg) per Q-particle capture. The temperatures relevant in helium- and carbon shell-burning are repre- 
sented by red bars and are marked with "He" and "C", respectively. 



reductions by up to an order of magnitude. The key rubidium 
and zirconium isotopes, for example, have undergone yield 
uncertainty improvements of roughly a factor of two. We have 
also found that s-process nucleosynthesis is more active when 
including the new 22 Ne+a reaction rates. While only small 
changes are found in the low-mass s-process path (A < 80), 
at higher masses production increases by up to a factor of 2. 
This is especially evident by counting the number of captures 



at the end of our network, yielding an increase of over 70%. 
Further calculations should be performed to study the effect of 
our new rates on lower mass AGB stars, while paying special 
attention to their effects on branchings in the s-process path. 

The Monte-Carlo method used in the present study to cal- 
culate the 22 Ne+a reaction rates has the distinct advantage of 
calculating the uncertainties in a robust and statistical mean- 
ingful manner. Although our rates include some of the sys- 
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7. (Colour online) Ratio of final abundances resulting from the new recommended 22 Ne(a,7) 26 Mg and 22 Ne(a,n) 25 Mg rates to those 



obtained from the old recommended rates. Points above unity (red line) represent a net increase in abundance, (a) At the end of core He- 
burning in a 25Mq star, the most significant abundances affected by the new rates are those of 26 Mg and the p-nuclei 74 Se, 78 Kr, and 84 Sr. 
(b) For AGB stars, higher mass nuclei are produced in larger quantities, as evidenced by the 'g' particle that monitors neutron captures above 
molybdenum. 



tematic uncertainties in the nuclear data, there are still open 
questions regarding the resonance properties that could affect 
the rates. Clearly, the remaining ambiguities in the nuclear 
data for the 22 Ne+a reaction rates need to be resolved. The 
discrepancies discussed here, by Koehler [60], and by Karakas 
et al. IU2I1 . make it difficult to assign some 26 Mg levels to 
22 Ne+a resonances. Furthermore, the E l f° = 831 keV reso- 
nance should be re-measured with high precision. More infor- 
mation should also be gathered on the structure of 26 Mg levels 
near the a-particle and neutron thresholds. Indirect methods 
such as particle transfer measurements are useful here, since 
the Coulomb barrier inhibits direct measurements. 
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